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We report the results of a first-principles study of dissociative electron attachment to H2O. The 
cross sections are obtained from nuclear dynamics calculations carried out in full dimensionality 
within the local complex potential model by using the multi-configuration time-dependent Hartree 
method. The calculations employ our previously obtained global, complex- valued, potential-energy 
surfaces for the three ( 2 Bi, 2 A\, and 2 B-z) electronic Feshbach resonances involved in this process. 
These three metastable states of H^O" undergo several degeneracies, and we incorporate both the 
Renner- Teller coupling between the 2 B\ and 2 A\ states as well as the conical intersection between 
the 2 A\ and 2 Bi states into our treatment. The nuclear dynamics are inherently multidimensional 
and involve branching between different final product arrangements as well as extensive excitation 
of the diatomic fragment. Our results successfully mirror the qualitative features of the major 
fragment channels observed, but are less successful in reproducing the available results for some of 
the minor channels. We comment on the applicability of the local complex potential model to such 
a complicated resonant system. 

PACS numbers: 34.80.Ht 




I. INTRODUCTION 



In the proceeding paper [l[ , referred to hereafter as pa- 
per I, we presented global representations of the three 
( 2 £>i, 2 Ai, and 2 i?2) complex- valued potential-energy 
surfaces of the metastable states of H2O - that under- 
lie dissociative electron attachment to water. This paper 
is concerned with the calculation of the cross sections for 
that physical process. Prior experimental and theoretical 
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have characterized the various breakup 
channels and determined the spatial symmetries of the 
three metastable electronic states of H2O - , the 2 E\, 2 A\, 
and 2 i?2 electronic Feshbach resonances, which arc re- 
sponsible for production of H~ and 0~. As explained in 
Ref. [IH and paper I, the energetically lowest H-f OER 
channel does not directly correlate with any of the three 
Feshbach states. We therefore conclude that OER pro- 
duction is due to nonadiabatic effects. 

We pursue this problem theoretically using a cou- 
pled Born-Oppcnheimer treatment of the nuclear mo- 
tion. The first task, which was described in paper I, 
is the construction of three-dimensional, complex-valued 
potential-energy surfaces for these three states, which 
have a negative imaginary component due to the finite 
probability of electron autodetachment back to H^O+e - . 
These complex-valued potential-energy surfaces, which 



are functions of the nuclear geometry q, are defined as 



V(q) = E R (q) - 



(i) 
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where En is the resonance position and T is the width 
of the resonance, which is related to the lifetime by r = 
1/r. (We use atomic units throughout this paper.) The 
present article, which we label paper II, is concerned with 
the use of these potential curves within the local complex 
potential (LCP) model HH, [H, [H, [H, H|] to calculate the 
nuclear dynamics leading to dissociation. The analysis of 
the dynamics yields the DEA cross section as a function 
of incident electron energy. 

We must account for two major nonadiabatic physical 
effects in calculating the quantum dynamics of the nu- 
clei. As described in paper I, the three potential-energy 
surfaces have several degeneracies that lead to coupling 
among them. First, the 2 B\ and 2 A\ states become mem- 
bers of a degenerate 2 II pair in linear geometry, and for 
this reason there will be Renner- Teller coupling between 
them. We expect this coupling to be relevant for DEA via 
the 2 A\ state, because the gradient of its potential-energy 
surface will cause the system to move toward linear ge- 
ometry after the electron attaches. Second, there is a 
conical intersection [r| between the 2 I?2 and 2 A\ states 
which leads to coupling between them. For this reason, 
as described in paper I, we constructed a set of diabatic 
2 i?2 and 2 A\ surfaces, along with a coupling term, which 
we use in the calculations presented in this paper. 

In Fig. [TJ we show the real parts En of the constructed 
potential-energy surfaces along a two-dimensional cut 
which includes the equilibrium geometry of the neutral 
(n = r 2 = 1.81a ; HO h = 104.5°). The degenera- 
cies that lead to the nonadiabatic effects listed above 
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FIG. 1: Real parts of resonance energies Er as constructed in 
paper I within Ci v geometry (ri = r?), plotted with respect 
to bending angle and symmetric stretch distance. 



can be seen in this figure. The two-dimensional cut 
depicted is that for which the two OH bond lengths 
are equal (r\ = ri), corresponding to Ci v symmetry. 
(In Ci v symmetry, the adiabatic and diabatic 2 A\ and 
2 i?2 surfaces coincide.) The backside of this cut lies 
at r\ = ?*2 = 1.81ao, which is the equilibrium value of 
the bond lengths in neutral H2O, and is marked with 
solid lines. The surfaces extend forward in Fig. [JJ along 
the symmetric stretch direction to geometries at which 
r\ = T2 = 2.7ao- The conical intersection comprises the 
set of points along which the 2 Ai and 2 B 2 surfaces in- 
tersect. The Renner- Teller degeneracy between the 2 B\ 
and 2 Ai states occurs at 0=180°. 

Although Fig. Q] shows only one cut of the potential- 
energy surfaces, and only their real part, it is useful for 
introducing certain features of these surfaces and the dy- 
namics that will result. Dissociative attachment via the 
lower 2 B\ and 2 A\ states leads primarily to the product 
H-+OH (X 2 n). The two OH bond lengths for such an 
arrangement are unequal, and therefore this product ar- 
rangement cannot be seen in Fig. [JJ However, we can 
see that at the equilibrium geometry of the neutral, the 
2 i?i surface is relatively flat with bend, while the 2 A\ 
surface slopes steeply downward toward linear geometry 
(0=180°). As a result, the dynamics beginning on the 
2 A\ surface will lead toward linear geometry, and we ex- 
pect that the Renner- Teller coupling between these two 
states will be more important for DEA via the 2 ^4i than 
via the 2 B\ state. 

The channel H2+0~ is the minor channel for DEA via 
the 2 -Bi and 2 A\ states, but the major channel for the 
2 i?2 state. We can see why this is the case from Fig. [1] 
the gradient of the 2 B 2 (2 2 A') surface leads downward 
toward the conical intersection, where the system may 
make a nonadiabatic transition to the lower surface and 
access the clearly visible H2+0 - well on the 1 2 A' (lower 
cone) surface. The 2 2 A' surface does not have a low 



energy asymptote in this geometry, instead correlating 

to 0-+H2 (ayj. 

The outline of this paper is as follows. In Sec. [Til wc 
summarize previous experimental and theoretical work 
on this problem. In Sec. IIII( we present the local complex 
potential model, which forms the foundation of our the- 
oretical implementation. The Hamiltonian for the rovi- 
brational nuclear motion of a triatomic molecule, and the 
additional terms which arise when the Renner- Teller ef- 
fect is included, arc described in Sec. IIVI In Sec. fVl we 
describe the multi-configuration timc-dcpcndcnt hartree 
(MCTDH) method, which we use to calculate the nuclear 
dynamics, and the formalism for calculating the DEA 
cross sections. In Sec. I VII we present the final results of 
this study: cross sections, as a function of incident elec- 
tron energy, resolved into the final rovibrational product 
states. 



II. PREVIOUS EXPERIMENTAL AND 
THEORETICAL RESULTS 

Dissociative electron attachment to water molecules 
has been the subject of previous experimental investi- 
gation, starting as early as 1930 0, and as recently as 
the past year (2006) [2(|. Early experiments on disso- 
ciative electron attachment to H2O focused mainly on 
the identification of the negative ion species formed, the 
measurement of the total cross sections, and the energy 
locations of the structures in the resonance process [y]. 
Buchel'nikova Q and Schultz p[ established that the 
main products of dissociative electron attachment to wa- 
ter are H~ and O - , with the production of 0~ being 
almost ten times smaller than that of H~ at lower ener- 
gies, but with 0~ dominating at higher electron-impact 
energies. 

Both Compton and Christophorou [f| and Melton Q 
carried out comprehensive studies of negative ion forma- 
tion in water and measured absolute cross sections for 
DEA. Three resonance peaks were observed. H~ produc- 
tion was observed at approximately 6.5 eV and 8.6 eV, 
with the second peak much less intense than the first. 
The species O - was observed in increasing intensities in 
three peaks at 7.0 eV, 9.0 eV, and 11.8 eV @. 

The species OH - is also observed in the dissociative 
electron attachment experiments, though at an inten- 
sity one order of magnitude below the minor 0~+H2 
channel, which is itself observed at an intensity approxi- 
mately one order of magnitude lower than the dominant 
H~+OH channel. Melton [|| argued that OH~+H was a 
true channel of dissociative electron attachment to H2O 
molecules, while in subseqent studies (e.g., Ref. [26|) it 
was argued that OH~ is produced by DEA to water clus- 
ters [H 2 0]„. The question of OH~ production has been 
reexamined in the recent experimental study of Fedor et 
a2.[2(|. These authors have concluded that, indeed, it is a 
direct product of dissociative electron attachment to wa- 
ter. This minor channel is not examined in the present 
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treatment, and no mechanism has, as yet, been advanced. 

The effects of isotopic substitution have also been an 
issue of some debate. The replacement of H2O by D2O as 
the molecular target has the effect of nearly doubling the 
reduced masses corresponding to OH (OD) bond motion. 
One would expect, at least in a simple one-dimensional 
picture, that the nuclear dynamics may be substantially 
altered by such replacement, and in particular, the time 
to dissociation is increased. A longer dissociation time 
allows a greater amount of electron autodetachment to 
take place; therefore, performing the same experiment 
with different isotopic variants provides information on 
the lifetime of the electronic state involved. The cross 
sections for DEA via both H 2 and D 2 were mea- 
sured, compared, and discussed in detail by Compton 
and Christophoreau Q ■ A smaller peak cross section for 
D~ production than for H~ production via the lowest- 
energy 2 B\ state was observed. On the basis of these 
results, these authors derived an approximate lifetime of 
2.1 x 10~ 14 seconds for the lowest-energy 2 B\ Feshbach 
resonance. We have already published an initial study 
of DEA via this state [l6|, which arrived at results and 
conclusions much different from those of Ref. p|. The 
calculated results yielded a higher peak cross section for 
D~ production via the 2 B\ resonance than for H~ pro- 
duction, and a similar energy-integrated cross section, in 
stark contrast to the results of Compton Q. The calcu- 
lations indicated a larger lifetime of 10.9 x 10~ 14 sec- 
onds for the 2 B\ state, and the nuclear dynamics that 
we calculated indicated that only a small portion of the 
dissociating anion flux is lost to autodetachment. 

The recent experimental results of Fedor et al. [2(| have 
substantially resolved this controversy. These authors 
obtain results different from those of Ref. reversing 
the trend in peak heights for H~ versus D - production 
via the 2 B\ resonance. They observe a higher peak for 
D~ production than for H~ production, which brings the 
current experimental and theoretical results into qualita- 
tive agreement. 

Although the peak heights provide considerable infor- 
mation about the physical process of dissociative elec- 
tronic attachment to water, further information is gained 
by resolving the angular dependence of the fragments 
produced, and the final (ro)vibrational state of the di- 
atomic fragment. A series of measurements by Trajmar 
and Hall [fl and Belie, Laudau, and Hall [9| revealed 
the energy and angular dependence of H~ in dissocia- 
tive electron attachment to H2O. The determination of 
the angular dependence aided the assignment of the spa- 
tial symmetries of the three resonant states, B±, A±, and 
B 2 , which had previously been misassigned. By resolving 
the kinetic energy of the H~ fragment, this experiment 
yielded information about the vibrational and rotational 
state distribution of the OH fragments. 

Curtis and Walker [lj| measured cross sections for dis- 
sociative electron attachment to D2O and obtained two 
important results. By measuring the kinetic energy of 
recoil of the D - fragments produced, these researchers 



established that both ground state OD ( 2 n) and excited 
OD ( 2 E) accompany the D~ anions produced within the 
third resonance peak, and that the three-body breakup 
channel D~+D+0 is observed toward the high-energy 
tail of the second peak. 

The experimental studies determined that there are 
three mctastable electronic resonance states of the H20~ 
anion, the 2 B\, 2 A\, and 2 i?2, which are primarily re- 
sponsible for dissociative electron attachment to water. 
These three electronic states correspond to the three 
peaks seen in the experimental cross sections. Although 
the third peak is not obvious in the H - cross sections, it 
is present, though much smaller than the first and second 
peaks. 

Several salient features of the early experiments sug- 
gest that the nuclear dynamics of this process may hold 
some surprises. For dissociative attachment through the 
2 Bi resonance, the cross section for producing H~ + OH 
is roughly 40 times larger at its peak than the cross 
section forproducing the energetically favored products, 
0~ + H2 The lowest-energy atom/diatom arrange- 

ment, H + OH~, is produced in even smaller quantities. 
In addition, the branching ratios for the different prod- 
uct states vary greatly depending on which Feshbach res- 
onance is formed by the attachment. These observations 
indicate that the products of this reaction are deter- 
mined by the dynamics of the process itself rather than 
by the energetics of the possible product channels, and 
that moreover those dynamics are different for each of 
the resonance states of the water anion. The detailed ex- 
periments of Belie, Landau and Hall Q in 1981 indicated 
that the dissociation dynamics involve correlated motion 
among multiple degrees of freedom. For instance, the 
channel producing H~ + OH through the 2 B\ resonance 
state is accompanied by extensive vibrational excitation 
of the OH fragment. 

Therefore, given the competition between dissociation 
channels and the observed product vibrational excitation, 
one expects that the dynamics of dissociative attachment 
to this molecule are intrinsically polyatomic, and can only 
be described theoretically by a treatment using the full 
dimensionality of nuclear motion. 

Compared with the large number of experimental 
measurements, detailed theoretical work on dissociative 
electron-water collisions has been relatively scarce. The 
paucity of theoretical work on DA stems from the fact 
that, in water, DA proceeds, not through tunneling shape 
resonances, but through Feshbach resonances that in- 
volve changes in the electronic structure of the target. 
Early theoretical work focused on the electronic struc- 
ture [ll[ and configuration-interaction [l2j calculations 
on various states of H20~ that are possible resonances. 
These calculations, together with experimental observa- 
tions, formed the basis of the assignment of the three Fes- 
hbach resonances that are responsible for electron-impact 
dissociation of water in the gas phase. 

Contemporary theoretical work has included ah ini- 
tio complex Kohn [l3| and i?-matrix [l4| calculations, 
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at the equilibrium nuclear geometry, of the resonances 
and excitation cross sections into low-lying dissociative 
electronic states. More recently, Gorfmkicl, Morgan, and 
Tennyson [l5| carried out i?-matrix calculations of disso- 
ciative excitation of water through the four lowest excited 
states (the 1,3 B± and 3,1 A\ states). A limited study of the 
effects of nuclear motion were included in that work by 
increasing one of the OH bonds while keeping the equi- 
librium HOH bond angle and the other OH bond length 
constant. The only theoretical work on the dynamical 
aspects of dissociative electron attachment to water are 
earlier classical trajectory analyses based on either repul- 
sive [27| or attractive [28| model resonace surfaces. 

We previously reported calculations of the cross sec- 
tions for dissociative attachment through the lowest- 
energy 2 By resonance [H, [l?} that incorporated a full 
quantum treatment of the nuclear motion of the resonant 
state. That study found good agreement with experiment 
for dissociative attachment through the lowest resonance 
state ( 2 Bi) of the water anion to produce H~, and it es- 
tablished that the associated dynamics are intrinsically 
polyatomic and thus cannot be described successfully by 
one-dimensional models. The present treatment super- 
sedes our earlier study and extends the treatment to in- 
clude the higher resonance states as well. 

We have recently presented a qualitative study of 
the potential-energy surfaces for the three Feshbach res- 
onances, which demonstrated that for these metastable, 
anion states, there exist numerous intersections and de- 
generacies within the adiabatic manifold. This study 
identified the conical intersection between the 2 A\ and 
2 i?2 states, as well as a novel degeneracy between the 2 Bi 
Feshbach resonance and a 2 £?2 shape resonance. This 
degeneracy defines a branch seam, and the two reso- 
nance energies are seen to comprise two components of 
a double-valued adiabatic potential-energy surface. This 
seam and the resulting dynamics may have an effect upon 
the three-body, H+H+0 _ cross section, although we do 
not include it in t he p resent treatment. Finally, in a sep- 
arate publication [l9( , we derived a "constant-eigenmode 
approximation" and used it to calculate the angular de- 
pendence of the H~ fragment production via the 2 B\ 
resonance. We found excellent agreement with the results 
of Belie, Landau and Hall @, and demonstrated that the 
observed angular dependence is a result of partial-wave 
mixing in the resonance-background coupling. 



III. LOCAL COMPLEX POTENTIAL MODEL 

We treat the nuclear dynamics of dissociative electron 
attachment within the local complex potential model. 
This model is concerned with the proper accounting for 
the decay of the resonant state, and its effect upon the 
nuclear dynamics. The LCP model includes the simplest 
such accounting, in which the decay rate is a local func- 
tion of the nuclear geometry. 



A. Feshbach partitioning and the nuclear wave 
equation 

The local complex potential model [2l|, [z| HI] , also 
known as the "Boomerang" model when applied to vi- 
brational excitation, describes resonance nuclear motion 
by an inhomogeneous Schrodinger equation and a com- 
plex, but purely local potential. It is perhaps easiest to 
derive by applying Feshbach partitioning [291 1 within the 
Born-Oppenheimer framework to derive a nuclear wave 
equation (24l.[25|. The derivation begins by defining a dis- 
crete (square-intcgrablc) approximation to the resonant 
electronic state, ipQ{f*e',Q) which depends parametrically 
on the nuclear coordinates q and which is unit-normalized 
with respect to integration over the electronic coordinates 
f e . One then defines the geometry-dependent Feshbach 
projection operator Q, which operates on the electronic 
degrees of freedom, 

Q(q) = \i>Q(q)]fyQ($)\, (2) 
and its complement P: 

P(q) = 1 - Q(q), (3) 

with P 2 = P, Q 2 = Q and PQ = QP = 0. (Brackets 
denote integration over the electronic degrees of freedom 
only.) Partitioning the full wave function for total energy 
E as \I> + = P^ + + , we can formally derive the 
following inhomogenous equation for Q^ + : 

(E-QHQ-QHP E _ p ^ p + i£ FHQ) Q ^ ^ 
= QHP*+, 

where H is the sum of the electronic Hamiltonian and 
nuclear kinetic energy, H = H e i + T^. 
In view of Eq. ([2]), we can write 

Q^ + {r- e -q) = i> Q (r e ;q) (5) 

The function £(g) describes the relative motion of the 
nuclei in the negative-ion resonance state. To derive an 
equation for t;(q), the first approximation that is made 
is the Born-Oppenheimer approximation: we neglect all 
non-adiabatic couplings arising from the operation of the 
nuclear kinetic energy upon the adiabatic basis. Then 
multiplying Eq. (0| from the left by il>Q(f e ;q) and inte- 
grating oven the electronic coordinates gives the nuclear 
wave equation, 

(E - V Q (q) - A(E) - T q -) = QH el PV+, (6) 

where 

V Q (q) = [^ Q \H el \^ Q } (7a) 

and 

A(E) ee QH el P — 1 PH e iQ. (7b) 

E - PH el P -Tq+ie 
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The real- valued potential Vq{q) is the expectation value 
of the electronic Hamiltonian with respect to the discrete 
state tpQ'i the additional, energy-dependent term A{E) is 
called the "level-shift operator" and is nonlocal in the 
nuclear degrees of freedom q, owing to the presence of 
the nuclear Green's function. The residue of this Green's 
function gives the level-shift operator A(E) a negative- 
definite imaginary component. 

In order to bring Eq. ([6]) into the form of the local com- 
plex potential model, it is necessary to make a local ap- 
proximation to the level-shift operator A(E), and also to 
approximate the driving term. The assumptions that un- 
derlie these approximations are well understood [3(1 HH . 
A local approximation to the level-shift operator yields 

V Q (q) + A(E)^E R (q)-i^f, (8) 

where En and T are the location and total width of the 
resonance. A first-order perturbation treatment (Fermi's 
golden rule) of the driving term yields [l9[ 



QH el P*+*^^ X »M = <l> V MU), (9) 

where Tq is the partial width for decay to the ground 
electronic state of the target, and Xv* IS the initial rovi- 
brational state of the target. 

The final working equation of the LCP model then 
reads 

(e - E R M + *M - Tf ) CM - ^x,M- 

(10) 

The location and widths of the various resonance states 
were obtained from configuration interaction and fixed- 
nuclei variational electron scattering calculations, respec- 
tively, as detailed in paper I. In the case of the 2 B\ reso- 
nance, which generally lies below its 3 i?i neutral parent, 
the resonance can only decay into the ground electronic 
channel. In that case, the total and partial widths, T and 
r , coincide and can be obtained by fitting the eigenphasc 
sum to a Brcit-Wigner form. For the higher resonances, 
a more elaborate fitting procedure is required to obtain 
the partial widths, as outlined in Ref. [19j and in paper 
I. 



B. Time-dependent formulation of the LCP model 

A direct solution of the differential equations of the lo- 
cal complex potential model can pose significant difficul- 
ties for problems with multiple degrees of freedom and, 
in such cases, a time-dependent formulation of the prob- 
lem can offer distinct computational advantages. Such a 
formulation can be made, as demonstrated by McCurdy 
and Turner [32j , by formally writing the solution of Eq. 
(HHD as 

^M = (E-H + ie)- 1 ^M0) (11) 



and writing the nuclear Green's function as the Fourier 
transform of the propagator for the time-dependent 
Schrodinger equation: 

/>oo 

CM = lim i / e^ E+i ^e- im <f> v M,0)dt 

~° ^ (12) 
= lim i / e l(E+le)t ^\q,t)dt, 
Jo 

where we define the time-dependent nuclear wave func- 
tion as 

^M,t) = e- mt ^M,0)- (13) 

The driving term <f> Vi (q, 0) of the LCP equation can 
thus be viewed as the initial value of a wave packet that 
subsequently evolves on the complex potential surface 
of the resonance anion. Since the potential surface is 
complex, the packet decays as a function of time until 
it effectively escapes the region of the surface where the 
width is nonzero. 



IV. TRIATOMIC JACOBI COORDINATE 
SYSTEM AND HAMILTONIAN 

The LCP model equations were solved in the coordi- 
nate systems depicted in Fig. O For the three internal 
degrees of freedom of this triatomic molecule, we employ 
Jacobi coordinate systems, which are depicted at the top 
of this figure. The Jacobi coordinate system on the left, 
marked "(a)," is used to analyze the OH+H arrange- 
ment; the one on the right, marked "(b)," is used for the 
H2+O arrangement. The vector r connects the nuclei of 
the diatomic. The vector R connects the center of mass 
of the diatomic to the third atom. R is the length of 
R, r is the length of r, and 7 is the angle between the 
R and r vectors. For (a), 7 = denotes a linear OHH 
configuration. 

In addition to the three internal degrees of freedom 
there are also the three Euler angles that orient the in- 
ternal or body-fixed (BF) frame with respect to the lab 
or space-fixed (SF) frame. The origin of both frames 
is the center of mass. The space-fixed Z-axis is always 
chosen to be parallel with the wavevector of the incident 
electron. For calculations with total rotational angular 
momentum J = 0, the Hamiltonian only operates on the 
internal degrees of freedom. For J / we must take the 
Euler angles into account, and we denote them by a, /3, £. 

The total (electronic plus nuclear, ignoring spin) an- 
gular momentum, J, and its projection upon the space- 
fixed Z axis, M, are quantum numbers conserved by the 
Hamiltonian. We also use the quantum number K to 
specify the projection of the angular momentum on a BF 
axis. K is not a conserved quantity and there is some 
flexibility in its definition. We use the "i?-embedding" 
scheme [H[ in which R is taken to be collinear with the 
BF Z' axis and the angular momentum number K is 
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FIG. 2: Jacobi coordinate systems used to analyze the 
OH+H, (a) and H2+O (b) arrangement channels and the "R- 
embedding" coordinate system with origin at the center of 
mass. Primed and unprimed axes refer to BF and SF frames, 
respectively. The BF X' Z' and X'Y' planes are both marked 
with a thin line circle and the SF XZ and XY planes are 
marked with dashed circles. The line of nodes is also drawn. 
The molecule resides in the BF X' Z' plane. 



quantized around this axis. With this convention, the 
Eulcr angles a and [3 arc the polar angles which orient 
the R vector with respect to the SF frame, and £ is the 
third Euler angle specifying orientation about the BF Z' 
axis. A schematic of the coordinate system is also shown 
in Fig. H 

We may write a general expression for the six- 
dimcnsional rovibrational wave function for a triatomic 
with specified J and M value as follows: 

Xui (R, r, 7, a, (3, Q = Y] D J MK (a, (3, Q Xu% r ' ^ , 
k Rr 

(14) 

where the basis of D J MK (a, (3, £) is the set of normalized 
Wigner rotation matrices (and BF angular momentum 
eigenstates) 



Dj IK (a,/3, C) 



2J+ 1 



D J MK (a,P,Q (15) 



such that 

da [ d(cos (3) [ d( D J MK (a, (3, QD J M t K ,(a : (3, () 
J-i Jo 

(16) 



and (HI 



1, which for the D J MK 



follow the conventions of 
is the same as that of 



In Eqs . (IFa 
Zhang 

Edmonds [35|. 

The standard [36j, [37| BF Hamiltonian for the ra- 
dial solutions °f tn i s expansion incorporates cou- 
pling among the different K values for a given total 
angular momentum J. The neglect of this coupling is 
termed the "coupled states" or "centrifugal sudden" (CS) 
approximation [38l [39j . and we employ this approxima- 
tion for our calculations, since the kinetic energy of the 
recoiling fragments is large compared to their centrifugal 
energy. The resulting Hamiltonian is thus diagonal in K 
and can be written 



K 2[i R dR 2 2fx r dr 2 



J' 



J' 



2jjL r r 2 2^ R R 2 



J(J+l)-2K 2 
+ „ '— + V(R, r, 7) 



J 



2 m R 2 

Id. d 

— sm7— 

sm 7 07 07 



(17) 



A - 



. 2 

sm 7 



where [i r and are the reduced masses in either de- 
gree of freedom and V is the (coupled set of) Born- 
Oppenheimer potential-energy surface(s) that we calcu- 
late. 



A. Inclusion of Renner-Teller coupling 

For dynamics beginning on the 2 Ai (1 2 A') resonance 
surface, the gradient of that surface will force the wave 
packet toward linear geometry, at which point this reso- 
nance state is degenerate with the 2 B\ resonance (see Fig 
P . The Rcnncr-Tcller effect jMESSllilElElliSIS 
will therefore couple these two components of the 2 I1 
state, and we modify the Hamiltonian of Eq. (fTT|) accord- 
ingly. 

The quantum numbers J(J + 1) and K in Eq. (|17[) are 
obtained as eigenvalues of the total angular momentum 
operators J 2 and J z >'- 

Jz> lD J MK {a,f3,Q Rr 

,K,t, \, ( 18 ) 



K[D J MK (a,f3,Q 



Rr 



etc., where J z > has a simple form in terms of derivative 
operators in (a, (3, £) [47|. Properly, the operators that 
appear in the Born-Oppcnhcimer Hamiltonian for the 
rovibrational motion of the nuclei, Eq. (|17[) . should be not 
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J 2 and J z i but the nuclear angular momentum operators 
R 2 and R z < , where 



Ji li 



(19) 



in which expression ij is an electronic angular momentum 
operator; the Hamiltonian [36| with this form is exact 
except for the omission of the mass-polarization term. 

The exact Hamiltonian [36[ introduces numerous new 
diagonal and off-diagonal (off-diagonal in K, electronic 
state, and both) coupling terms to the triatomic Hamil- 
tonian. The term that is most commonly labeled 
the Renner- Teller coupling comes from the j 2 term in 



1 1 

+ 



2ji r r 2 2fi R R 2 J sin 7 



K- 



1 



1 



2fi r r 2 



Rl 



2(i R R 2 ) sin 2 7 



2\x r r 2 2\i R R 2 



K 2 - 2Kf z , + I 2 



■ 2 

sm 7 



(20) 



It is the 2Kl z i term which couples the two components 
(sine and cosine, 2 B\ and 1 2 A' ) of the 2 H state at 
linear geometry. At such geometries the operator l z i is 
diagonalized by 



iy (lj)A' ± ilpBl) = ±1 X (lj) A , ± ilpBl) 



(21) 



The matrix elements of l 7 j may either be computed |4ll . 
5H, l46| . or approximated by their values at linear geom- 
etrv [42l. |43| . We take the latter route, i.e., we assume 
that Eq. (|21|) holds everywhere. This approximation has 
little effect on the dynamics because only near linear ge- 
ometry docs the coupling become large. We perform our 
Renner-Teller calculations in the (l z > = ±1) diabatic ba- 
sis because it allows us to incorporate the boundary con- 
dition in 7 using the "JT-legendre" discrete variable rep- 
resentation [37| . With this assumption, for a given value 
of K, the (l z i = ±1) diabatic states have R z / = K ± 1. 
The kinetic energy operator in Eq. (f^Uj) is diagonal in this 
diabatic basis. The coupling then arises from the elec- 
tronic Hamiltonian, which is not diagonal in this basis. 
The electronic Hamiltonian in this basis takes the form 



V = - ( Va ' + Vb i Va ' ~ Vb > 
2\Va>- V Bl Va- + V Bl 



(22) 



i.e., the diabatic states are degenerate. When K = 0, 
there is no Renner-Teller effect, since the coupling term 
in Eq |2"01 vanishes. 



V. THE MULTICONFIGURATION 
TIME-DEPENDENT HARTREE METHOD 

The Multiconfiguration Time-Dependent Hartree or 
MCTDH [H, [H, \m HI method is an efficient adaptive 



scheme for propagating quantum-mechanical wave pack- 
ets for systems with multiple degrees of freedom. We use 
this method to perform the propagation in Eq. (fTU|) . We 
use the implementation within the MCTDH package [52l|. 
a freely available suite of codes built at the University of 
Heidelberg, Germany. 

In the MCTDH method, as in other methods developed 
for solving the time-dependent Schrodinger equation, we 
start with a time-independent orthonormal product basis 
set, 



{X^(<li)-X^(qf)}, J,. 1---A, 



(23) 



for a problem with / degrees of freedom and nuclear co- 
ordinates labeled qi,...qf. For computational efficiency, 
the basis functions Xj are chosen as the basis functions 
of a discrete variable representation (DVR) 

The central idea of the MCTDH technique is the repre- 
sentation of the nuclear wave packet as a sum of separable 
terms, 

^Mt) = E - E A n-jf(t) II Pi? («»>*), ( 24 ) 

with n K -c N K . Each "single particle function" (or SPF) 

ifjj(q K ,t) is itself represented in terms of the primitive 
basis, 



<P 



£ ) fe,*) = ECi(%rfe) 



(25) 



One can determine equations of motion for the param- 
eters c 2 - K ^ (t) and Aj 1 _,,j f (t). Since both the coefficients 
and the single-particle functions are timc- 



-4 



n—jf 

dependent, the wave function representation is made 
unique by imposing additional constraints on the single- 
particle functions which keep them orthonormal for all 
times (Hof . 

The evaluation of the Hamiltonian matrix, which 
must be carried out at every time step, may be expe- 
dited (49l [50l ] if the Hamiltonian can be written as a sum 
of products of single-coordinate operators. The MCTDH 
package [52( includes a utility which performs a fit of a 
given potential to a separable representation of this form. 
Details can be found in Beck et a/.|50j]. All potential- 
energy surfaces used in the current calculation were rep- 
resented in this manner, using this utility to fit them 
specifically for each choice of the DVR grids. 

For calculations on the electronically coupled 2 A\ and 
2 i?2 states, the underlying DVR is the same for each elec- 
tronic state, but each electronic state has its own set of 
single-particle functions <p) ■ This is referred to as the 
"multi-set" formalism, as opposed to "single-set." The 
Renner-Teller coupled 1 2 A' - 2 B\ calculation is per- 
formed under the single-set formalism. 
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A. Complex absorbing potentials 

The sine DVR bases in the r and R degrees of free- 
dom incorporate standing wave boundary conditions at 
their edges. Therefore, when the dissociating wave packet 
reaches the end of the DVR grid, it must be absorbed 
to prevent unphysical backward reflections. To this end 
we include an artificial negative imaginary component 
to the surface called a "complex absorbing potential" or 
CAP S [H: 



V ( 



CAP = 



iij(R 







(R < R c ) 
(R > R c ); 



(26) 



a similar expression for the CAP in the r degree of free- 
dom also applies. Formally, the CAP's provide the +ie 
limit in Eq. (fT2l) . 

We use a value for 77 equal to 0.007 hartree, and place 
R c three bohr before the end of our grid, except for the 
1 2 A' calculations for H2 and D2, for which we use a 
strength of 0.0018 hartree and a value of R c five bohr 
before the end of the grid. 



B. Dissociative attachment cross sections from 
outgoing projected flux 

The cross sections for dissociative attachment can be 
calculated directly from the time-propagated wave packet 
by computing the energy-resolved, outgoing projected 
flux. The energy resolution is achieved by Fourier trans- 
form and a final state resolution is achieved by the intro- 
duction of appropriate projection operators. For DEA to 
a specific final rovibrational state labeled by rotational 
(j) and vibrational {v) indices, we use the projection op- 
erator 



r 



Xju 
r 



(27) 



The flux operator, which measures the flux passing 
through a surface defined by R = R c , is defined as 



F = i[H,h(R-Rc)}, 



(28) 



where h is a hcavisidc function. The energy-resolved pro- 
jected flux is then given by 



— / dt d£ 
2tt J Jo 



b Vi \e^ H - E ^P jv FP jv e-^ H - E ^'\4> Vi ). 



(29) 



The MCTDH package [52| includes a utility which com- 
putes the outgoing projected flux. In the actual calcula- 
tions, the flux operator appearing in the equation above is 
replaced by an expression involving the complex absorb- 
ing potential, Eq. ([26|) . This formulation of the flux oper- 
ator is very convenient numerically and entirely equiva- 
lent to the traditional formal definition of the operator in 



this context, in the limit that the CAP does not perturb 
the propagating wave packet beyond first order, which 
in the present case, given the nuclear masses, holds as a 
good approximation. For more details on this CAP flux 
formalism see Refs. [Hoi. IHll l56j . 

The resulting energy-resolved projected flux is that as- 
sociated with the time- independent solution of the driven 
Schrodingcr equation of the LCP model in Eq. HI 



In terms of Fj v , the DEA cross section is [ij 



DEA 



4tt 3 



2 



(30) 



(31) 



For the H~+OH channel, an additional factor of two is 
included in Eq. (f3Tj) to account for the fact that in a 
given calculation we perform the flux analysis for only 
one of the two H~+OH arrangements, namely the one 
for which the Jacobi coordinates are appropriate. 

The definition of the rovibrational states Xjv is com- 
plicated by the ion-dipole interaction of the fragments. 
In our earlier study DEA to water via the 2 B\ Fesh- 
bach resonance [17], we attempted a complete final state 
analysis, and projected upon pendular (restricted rotor) 
states [53], not free rotational states, and assumed that 
these pendular states evolve adiabatically to their free 
rotational state asymptotes. This analysis did not yield 
any major insight, and so for the present calculations we 
simply project upon free rotational states. As a conse- 
quence, there is a small error in our final state resolution, 
but the magnitude of this error will span a range of states 
and range of energies approximately equal to the magni- 
tude of the ion-dipole interaction at the edge of our grid, 
which is small compared to the kinetic energy spread of 
the fragments. 



C. The DVR bases and other MCTDH parameters 

In most of the calculations reported here, we used 
DVR primitive basis sets for all internal degrees of free- 
dom, choosing the standard sine DVR[f3(| for the r and 
R degrees of freedom and, for J = 0, the Legcndrc 
DVRj50| for 7. For J > 0, as previously discussed, the 
DVR for 7 must be modified to account for singular- 
ities in the Hamiltonian [see Eq. pT)) ] due to the term 
K 2 /sv^_(jY This is done by using an extended Legendre 
DVR [13, HI, [1^] , which is implemented in the Heidelberg 
MCTDH package (H. 



D. Initial states 

The initial rovibrational states \vi °f Eq. (|9|) w ere ob- 
tained via relaxation and improved relaxation [5l| as im- 
plemented within the MCTDH package In relax- 
ation runs, an initial guess Xg(<?:0) for the ground state 
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is propagated in imaginary time, which yields the ground 
state Xo(o)- 

Xg(q,T)=e- Hr Xa (q,0) — ► Xo(q) (32) 

r — >oo 

In improved relaxation runs, the propagation of the SPF 
expansion coefficients c|*]- k (t) of Eq. |25|) is performed via 
Eq. (|52")) . but the configuration coefficients A^. (t) are 
obtained anew at each time-step via a Davidson diago- 
nalization. 

For the two-state 2 B2- 2 A\ calculations, the wavefunc- 
tion is represented in the diabatic basis, each compo- 
nent of which has an expansion of the form of Eq. ([M|) . 
and different sets of time-dependent single-particle func- 
tions <pW (q K ,t). Since the adiabatic-to-diabatic transfor- 
mation angle is not constant with geometry, the single- 
particle functions which represent the initial state will be 
different in the diabatic basis than in the adiabatic basis 
(in which they would be identical to the single-particle 
functions of the relaxation run). For this reason, an it- 
erative technique (50l | is employed to minimize the er- 
ror between the diabatic representation and its adiabatic 
representation. 



VI. CALCULATED CROSS SECTIONS FOR 
DISSOCIATIVE ELECTRON ATTACHMENT TO 
WATER 

We will present cross sections for dissociative electron 
attachment to water into the three different atom-diatom 
arrangements that are present as asymptotes of the 2 B\ 1 
2 Ai, and 2 B 2 Feshbach resonances, 

{H 2 + Cr 3.56 cV 

H + OH (X 2 U) 4.35 eV 
H + OH* ( 2 S) 8.38 eV, 

resolved into the final rovibrational states of the diatomic 
fragment, as a function of incident electron energy. The 
final state resolution allows us to determine the kinetic 
energy of the diatomic fragment. Therefore, we are able 
to calculate cross sections as a function of both the inci- 
dent electron energy and the kinetic energy of the re- 
coil, which data we may easily compare with experi- 
ment. These two dimensional data provide a comprehen- 
sive view of dissociative attachment via each of the reso- 
nances. We calculate the degree of rotational and vibra- 
tional excitation, and show how these quantities change 
with the incident electron energies. 

We have obtained converged cross sections for all chan- 
nels considered, with two exceptions. For the Renner- 
Teller coupled 2 A\ (I 2 A') and 2 B\ states, we have been 
unable to obtain a non-zero result for the minor H2+O - 
channel, as well as its deuterated counterpart. For the 
production of H~ from the 2 B 2 state, coupled to the 2 A\ 
state via the conical intersection, our calculations are not 



fully converged, although we do obtain total cross sec- 
tions and, for the OH ( 2 £) fragment, final state resolu- 
tion. 

In general, our results arc in qualitative, though not 
quantitative, agreement with the experimentally mea- 
sured cross sections for the major product arrangments 
observed in DEA via each resonance state: H~+OH 
from the 2 B\ and 2 A\ resonances, and H2+0~ from the 
highest-energy 2 B 2 resonance. We have had difficulty 
obtaining results for the minor channels, failing to repro- 
duce the experimental result for the magnitude of the 
cross section for production of H2+O - via the first two 
resonances, and not being able to fully converge the cal- 
culation for production of H~+OH ( 2 n and 2 E) via the 
2 Bi resonance. 

Total cross sections calculated for the two anion- 
diatom arrangements, and for both D2O and H2O, are 
presented in Fig. [31 along with recent experimental re- 
sults from Fedor et al.^f^. Calculated and experimental 
peak heights and locations arc collected in Table fl] In 
Fig. [31 the experimental data are internormalizcd but not 
absolutely normalized; therefore, for the purposes of com- 
parison, we normalize the experimental peak for the pro- 
duction of H~+OH via the lowest-energy 2 B\ resonance 
to Compton and Christophoreau'sQ result of 6.5xl0~ 19 
cm 2 . The calculated curves are obtained by summing 
the cross sections into the individual rovibrational states 
Xjv of the ion + diatom arrangement. Therefore, the 
three-body cross sections are neglected. The recent ex- 
perimental results of Ref. [2(| do not resolve the kinetic 
energy of the atom-diatom recoil and therefore do not 
distinguish between the two- and three-body DEA cross 
section. Thus, to the degree that three body breakup is 
important, our calculated results cannot be compared di- 
rectly with these experiments, for incident electron ener- 
gies which exceed the three-body dissociation thresholds 
of cither 8.04cV (H+H+0~) or 8.75cV (H+H~+0). 

We can draw the following conclusions from Fig. [3] and 
Table [H First, it is clear that the entrance amplitude for 
the 2 .E?i resonance has been overestimated by our present 
study, because the magnitude of the DEA cross section 
via this resonance is entirely controlled by its entrance 
amplitude, and we have overestimated the experimental 
peak height by nearly 60%. Also, similar to the result 
of our previous studv|16l.[T7|. we overestimate the energy 
at which the H~ cross section via the 2 B\ state peaks 
by about 0.4eV. The calculated peak location, 6.87eV, 
is larger than the vertical transition energy for the 2 B\ 
resonance as defined by our potential-energy surface con- 
structed in paper I, which is 6.63eV. This value is ob- 
tained through a configuration-interaction treatment of 
the resonance; using complex Kohn scattering calcula- 
tions, we obtained a value of 6.09eV. The comparison 
between the calculated and experimental peak locations 
indicates that the physical value of the vertical transition 
energy for the 2 i?i Feshbach resonance is probably about 
6.2eV, nearer to the complex Kohn result. 

Similar observations apply to the comparison between 
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Incident electron energy (eV) Incident electron energy (eV) 




Incident electron energy (eV) Incident electron energy (eV) 



FIG. 3: Cross sections for production of HP or D~ (left) and (right) as a function of incident electron energy, summed 
from different MCTDH calculations. Top, cross sections from H2O; bottom, cross sections from D2O. Experimental results of 
Fedor et al. [2(| included for comparison. The experimental data, which do not have absolute normalization, are normalized to 
agree with Compton and Christophoreau's [f| H~+OH peak height for DEA via the 2 Bi resonance at 6.5 x 10 cm . 



the calculated and experimental results for DEA via the 
2 A\ resonance to produce H~+OH; the calculated peak 
height is too large and located at a higher energy than 
is experimentally observed. Therefore, it is possible that 
we have overestimated the entrance amplitude and the 
vertical transition energy for this resonance as well. The 
transition energy as defined by our potential-energy sur- 
face is 9.01eV; from complex Kohn calculations we ob- 
tained a lower result of 8.41cV, which is probably closer 
to the physical value. However, as we explain further 
below, the disagreement in magnitude and location be- 
tween the calculated and observed results may indicate a 
breakdown of the local complex model for DEA via the 
2 Ai resonance. 

The data in Fig. [3] and Table |T] indicate that, 
while the calculations overestimate the cross sections for 
H~+OH production, they evidently underestimate those 
for H2+0~ production. H2+0~ is the major channel for 
DEA via the 2 £?2 resonance. As explained in Ref. [l|[, the 
presence of this channel is entirely due to nonadiabatic 
coupling between the upper 2 Bi (2 2 A 1 ) resonance and 
the lower 2 A 1 (1 2 A') resonance via their conical intersec- 



tion. As we will show, the magnitude of this cross section 
is determined by active competition between different 
product arrangements, the dynamic effects of both the 
real and imaginary components of the surface, as well as 
the conical intersection dynamics. We regard the agree- 
ment with experiment that we have obtained to be quite 
good, considering the complexity of the system. We note 
that the location of the calculated peak for 0~ produc- 
tion which we have calculated (11.75eV) agrees well with 
the experimental value (11.8eV). The peak location may 
be contrasted with the vertical excitation energy which 
was calculated in paper I to be 12.83eV. The peak max- 
imum is a full 1 eV below the vertical transition energy, 
which difference reflects the influence of autodetachmcnt 
upon the nuclear dynamics. The large autodetachment 
probability weights those components of the propagating 
wave packet which are closer to the product arrangement, 
i.e., lower on the potential-energy surface, and results in 
a breakdown of the multidimensional reflection principle. 

The product channel H2+0~ is the minor one for DEA 
via the first two resonances, 2 B\ and 2 A\. We have failed 
to reproduce the corresponding experimental results; our 
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TABLE I: Peak cross sections (a, in units of 1CP 19 cm 2 ) and 
peak locations (E, in eV) calculated for DEA via the three 
resonances, . Experimental peak locations are taken from the 
data of Fedor et al, except where noted. 





i, n 1 r*i list pH 

2 Bi 2 Ai 2 B 2 


FiYTipri in Ptit 

JJAJJCl llilCULi 

2 Bi 2 A x 2 B 2 


OH + BT 

(x 2 n) 


a 103.7 41.4 2.61" 
E 6.87 8.74 11.54 


65" 13" 
6.4 8.4 


OH + H~ 

( 2 E) 


a 3.67" 
E 12.68 




H~ 
(total, 2 B 2 ) 


a 5.21" 
E 12.61 


? 

11.8 C 


H 2 + 0~ 


a 0.121 1.87 
E 7.62 11.75 


1.3 d 3.2 d 5.7 d 
7.1 9.0 11.8 


OD + D" 

(x 2 n) 


a 124.4 41.6 1.45° 
E 6.93 8.87 11.93 


52" 6" 
6.4 8.5 


OD + D" 

( 2 E) 


a 1.96" 
E 12.86 




D" 
(total, 2 B 2 ) 


a 2.60" 
E 12.57 


? 


D 2 + O" 


a 0.0242 1.97 
E 7.63 12.10 


? ? ? 
7.1 9.0 12.0 



"Calculation not converged for H +OH arrangement via B 2 res- 



onance 
b 



Compton and Christophoreau, Ref. 0| 
c Jungcn, Ref. Il2| | 
d Mclton, Ref. W 



calculations produce a very small cross section for DEA 
via the 2 B\ resonance to produce this channel, and zero 
cross section for 2 A\. We regard nonlocal effects to be a 
prime candidate for the physical origin of this channel for 
the 2 Bi state; for the 2 A\ state, we suspect that three- 
body dissociation into H+H+O", which is not treated 
in the present study, may play a significant role in this 
channel. 

The probability of a dissociative attachment event is 
neatly divided into distinct probabilities for attachment 
and survival by the local complex potential model. The 
norm of the driving term in the driven Schrodingcr equa- 
tion of the LCP model, Eq. (flO|) , corresponds to the prob- 
ability per for electron attachment, weighted by the en- 
velope of the initial vibrational state. We define this 
quantity as the attachment width Ta, 



r_ 4 = 2n {(j> Vi \<j> Vi ) , 



(33) 



and list the values of Ta for each of the resonances in 
Table HI 

Once the electron has attached, the loss of flux via 
the imaginary component of the complex- valued surfaces 
determines the survival probability of the anion state. 
The survival probability may be calculated by integrating 
the flux Fjv(E) over energy: 



J dE F AE) 



(34) 



The calculated survival probabilities are also listed in Ta- 
ble |TTJ The survival probability for the lowest 2 B\ res- 



TABLE II: Attachment widths and survival probabilities cal- 
culated for the three resonances using Eqs. (|33l) - (I35I) . 



Attachment width 



( xlQ- 



Survival 
probability 





2 B X 2 A x 2 B 2 


2 B X Mi 2 B 2 


H 2 
D 2 


3.47 4.33 4.74 
3.54 4.19 4.84 


93.8% 65.1% 21.5% 
91.6% 57.2% 13.1% 



onance state is near 1, and therefore for this resonance 
the magnitude of the cross section is controlled by the 
attachment probability. For DEA via the 2 A\ state, the 
cross section is lowered by the effect of autodetachmcnt, 
though once attached the electron is more likely to sur- 
vive to dissociation. For DEA via the upper 2 Bi reso- 
nance, the large majority of the attached wave packet is 
lost to autodetachment; therefore, the variation of the 
lifetime of this state with nuclear geometry plays a large 
role in the dynamics. 

We calculate the average degree of rotational and vi- 
brational excitation, as well as the average kinetic en- 
ergy of the anion recoil, for each of the final channels 
by weighting the survival probability by the quantity of 
interest, 



(Ekin) = 



(f) 



y. 



M,, 



_ E^vfdE F 3V {E) 

Psurv {.^fi 1 4*vi ) 

V „.,,,/ • /■;,.:/■■ 

Psurv (fivi | ) 

■{Einc — Ej V ) J dE Fj V {E) 



(35) 



In the third line of Eq. (|35|) , £!j nc is the incident electron 
energy, Ej V is the energy of the final state relative to the 
ground vibrational state of H2O, and M to tai and Mdiatom 
are the masses of the original triatom and the diatomic 
fragment, respectively, so that the resulting quantity is 
the kinetic energy of the anion recoil in the laboratory 
frame. We present these results in Table IIII1 along with 
experimental data on the anion recoil kinetic energy from 
Fedor et al. [Ij. 

Our calculated average values of the anion recoil ki- 
netic energy, agree to varying degrees with the 
results of Ref. [2(| ■ Our calculated values for average ki- 
netic energy release for the production of H~ from the 
2 B\ and 2 A\ resonances are much larger than observed 
by these authors, but closer to the values measured by 
Belie et al. @. As described in paper I, the potential- 
energy surfaces which we have constructed for these res- 
onance states reproduce the energetics of the two-body 
asymptotes very well (to within 0.08eV for the ground vi- 
brational state). Errors in the present results may there- 
fore only come from errors in the vertical transition en- 
ergies, or a misrepresentation of the dynamics prior to 
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TABLE III: Expectation values of final vibrational quantum 
number, {u), and that of angular momentum quantum num- 
ber, (j 2 ), of diatomic fragment, as well as the expectation 
value of the kinetic energy of the anion recoil, {E kin ), for 
each resonance, as calculated with Eq. (|35|l . Average kinetic 
energy determined by the experimental method of Ref. [20l |. 
final column. 



Diatomic fragment 


(u) {f) (E kln ) (E kln ), exp't." 


' 2 B 1 

tr+OH 

H 2 + 0- 
D" + OD 
D 2 + 0- 


1.28 107 2.04eV 0.96eV 
(1.5eV) 6 
3.52 412 0.154 0.12 
1.75 240 1.89 0.70 
3.34 821 0.125 0.14 


(1 'A') 
H" + OH (A'=0) 

D" + OD (A=0) 


2.11 121 3.35 1.55 

(2.5 eV) 6 
2.98 221 3.16 1.20 


Z B 2 (2 M') 
ET + OH ( 2 n) 
H" + OH ( 2 £) 

H 2 + 0- 
D" + OD ( 2 n) 
D - + OD ( 2 £) 

D 2 + 0- 


C C C 

4.69 c 439 c 3.45 c 

7.75 405 0.413 0.57 

c c c 

6.63 c 819 c 3.35 c 
13.0 725 0.684 0.79 


Data from Ref. l20t . except where noted. 



c The calculation for H production via E>2 is not converged. 



breakup. As discussed above, it is likely that our vertical 
transition energies for the 2 B\ and 2 A\ resonances are 
too high, perhaps by as much as 0.4eV relative to the 
proper physical values. Most of this excess energy may 
be transmitted to the kinetic energy of the H~ recoil, 
due to the small mass of hydrogen relative to the H 2 
molecule. While our calculated results exceed the exper- 
imental result of Ref. [2(| by more than 0.4eV, they are 
within 0.4eV of the Belie et al. value for the 2 Bi state 
and 0.85 eV for the 2 A\ state. Fedor et al. comment that 
the values obtained for the kinetic energy release of H~ 
via the 2 B\ and 2 A\ resonances by Belie et al. ^ "may 
be considered as more accurate." The discrepancy with 
Belie et al. for 2 B\ result supports our recommendation 
that the physical transition energy for the 2 B\ state be 
taken to be approximately 6.2eV. The maximum kinetic 
energy release at our calculated peak locations of 6.87 
and 8.74eV for the 2 B\ and 2 A\ states, is, respectively, 
2.38 and 4.15eV. Our results are therefore very near the 
maximum values and reflect the small average degree of 
vibrational excitation which we calculate. 

Our values for the average kinetic energy release of the 
major 0~ fragment from the 2 i?2 resonance agree much 
better with the results of Fedor et al.; we again underesti- 
mate the experimental result, but only by 28% and 13%, 
respectively, for the nondeuterated and deuterated tar- 
get. For this channel, the degree of excitation of the H 2 
(D 2 ) fragment is large, and therefore the kinetic energy 
of the atom-diatom recoil is less than its maximum al- 



lowed value. The maximum kinetic energy release at our 
calculated peak (11.75eV) is 0.91eV, more than twice our 
result for the average value. Therefore, more energy goes 
into the rovibrational excitation of the H2 fragment than 
into the kinetic energy of the recoil. 

We do not calculate the three-body dissociative elec- 
tron attachment cross section, i.e., 



H 2 



H + H + O- 
H~ + H + O 



8.04 eV 
8.75 eV. 



(36) 



The complex absorbing potential flux formalism [50 . l5ll . 
[56| which is employed within the MCTDH implementa- 
tion (52| is not appropriate for the three-body breakup 
channel, at least when used in conjunction with the Ja- 
cobi coordinate systems used here. We do, however, pro- 
duce rigorous results for the two-body channels by pro- 
jecting upon the bound rovibrational final states as in 
Eq. (f2"T|) and summing. 

Our surfaces, as described in paper I, are not designed 
to reproduce the dynamics leading to three-body disso- 
ciation either. Due to our neglect of the shape/Feshbach 
resonance intersection on the 2 _B 2 manifold, which is 
a true characteristic of the physical system and which 
leads to the double- valuedness of the physical surface, 
we cannot accurately represent the dynamics leading to 
the three-body dissociation channels with our single 2 _B 2 
surface. The 2 £? 2 manifold is coupled to the 2 A\ state 
in the three-body region by the conical intersection, and 
therefore it is possible that this omission affects the dy- 
namics via the 2 A\ state as well. 

It is clear that we may only rigorously compare our 
results with experiment for the two-body channels. The 
comparison is complicated by the fact that the exper- 
imental results are sometimes not final-state resolved, 
such as those presented in Fig. [31 and in any case al- 
ways incorporate a finite resolution in determining the 
energy of the incident electron energy and the kinetic en- 
ergy of the recoil. The energetics of the asymptotes of 
the physical surfaces, which are mirrored very well by our 
constructed potential-energy curves, dictate that for the 
lowest 2 B\ state the three-body channels are closed, but 
that for DEA via the other two resonances, at least one 
three-body channel is open. 

In the following subsections, we give the principal find- 
ings of the nuclear dynamics calculations for each chan- 
nel that was studied. Further details are given in the 
EPAPS archive[6(|- Most calculations were carried out 
for the ground vibrational state and for total rotational 
angular momentum J = (or, in the case of the Renner- 
Teller coupled 1 2 A'- 2 Bi calculations, R z > = 0). Rovi- 
brationally excited initial states were also examined, and 
these are listed in the descriptions of the individual cal- 
culations which follow. 
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A. Dissociative electron attachment via the 2 B 1 
state 

We have previously [H, [l7| performed a calculation on 
the 2 B\ state, which is superseded by the present treat- 
ment. We perform the calculation using the one uncou- 
pled 2 B\ potential-energy surface. We have confirmed 
that Rcnner- Teller coupling to the 2 A\ state at linear ge- 
ometry has a negligible effect on the dynamics, at least 
for DEA via the ground rovibrational state of the target. 

The treatment in our previous study [l6l . \v\ was not 
able to resolve the cross section in the minor 0~+H2 
channel (A of the major channel), due to deficiencies 
in the potential-energy surface. While the present study 
does obtain converged _ +rl2 cross sections, they are 
two orders of magnitude below the observed cross sec- 
tions, and peak at energies well above the experimental 
peaks; it is therefore clear that we have not represented 
the dynamics leading to this minor channel accurately. 
The failure in this regard may be due to small deficien- 
cies in the potential-energy surface, or to the presence of 
significant nonlocal effects in this minor channel. 

In the previous study, we reproduced the magnitude of 
the OH+H~ cross section to within a few percent. We 
continue to regard that close agreement to be essentially 
fortuitous. Apart from an overall scale factor, the present 
calculations for production of H~ reproduce the shape 
and energy dependence of the experimental results very 
well. It appears that the current value calculated for the 
width of the 2 B 1 resonance at the equilibrium geometry 
of the neutral, 10.31mcV, is larger the physical value, the 
latter being closer to our previously calculated value of 
6.0mcV. It is likely that a more accurate description of 
electron correlation than we could include in the complex 
Kohn calculations of paper I is required to reproduce the 
resonance wavefunction of the 2 B\ state. 



Production of OH (X 2 U) + H~ 



via Bi state 



This is the dominant channel for DEA to H2O, hav- 
ing a peak cross section of approximately 6xl0~ 18 cm 2 . 
Cross sections as a function of incident electron energy 
are shown in Fig. [U We calculate a peak cross section 
of 10.35 x 1(T 18 cm 2 at 6.87eV. The magnitude of this 
cross section is larger than the experimental value (6.6 
x 10~ 18 cm 2 ), and the location of the peak is displaced 
upward by 0.4eV. For this resonance, autodetachment is 
nearly negligible. Therefore, the excess in the magnitude 
which we calculated reflects the fact that the calculated 
width values, and hence the entrance amplitudes, are too 
large. 

The cross sections in Fig. 0] are very similar in shape 
to those produced previously [H, QjJ, though they are 
larger in magnitude. At low incident electron energies, 
the first vibrational state is produced exclusively, and 
subsequent vibrational states have sharp onsets. As the 
degree of vibrational excitation increases beyond the first 
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FIG. 4: Cross sections for production of H _ +OH (y), top, 
or D~+OD (u), bottom, from 2 Bi state as a function of in- 
cident electron energy. Total, thick line; vibrational states 
(ground) through 5, thinner lines. 



few quanta, the magnitude of the cross section decreases. 
Despite the fact that the first five excited vibrational 
states are clearly visible in Fig. 0] the average number of 
quanta excited is only 1.28. The average kinetic energy 
release, therefore, is near its maximum value of 2.38eV. 
The degree of rotational excitation calculated for this 
state is relatively low ((j 2 )=107 for H2O). 

A two-dimensional view of the data is provided in 
Fig. where the kinetic energy of the anion recoil, as de- 
termined by a full final-state resolution of the products, is 
plotted versus incident electron energy; the contour lines 
indicate the magnitude of the cross section. The kinetic 
energy of the anion recoil in the laboratory frame is 
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total 



(37) 



as in Eq.j35|). This figure shows that the degree of ro- 
tational excitation for production of both H~+OH and 
D~+OD is small compared to the vibrational spacing of 
the OH fragment, because there are separate lobes cor- 
responding to each vibrational state. For the deuterated 
case, the lobes are thicker and closer together. The thick 
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FIG. 5: Top, cross section for production of H~+OH from 
2 B\ resonance as a function of incident electron energy and 
H~ fragment kinetic energy, unshifted, with the physical value 
of the maximum kinetic energy available plotted as bold line; 
bottom, deuterated. The physical value for the maximum 
kinetic energy is slightly lower than the value corresponding 
to our calculated surfaces. Contours every 2 x 10~ 17 cm 2 
eV" 1 . 
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FIG. 6: Top, production of H +OH via the 2 Bi resonance 
at different incident electron energies, as a function of H~ 
fragment kinetic energy, unshifted, on arbitrary and different 
scales for each incident electron energy. Bottom, D~+OD. 
Calculated results have been broadened using a 150meV 
linewidth, consistent with the plotted experimental results 
from Belie, Landau and Hall 0]. 



line in this figure corresponds to the maximum kinetic 
energy available, as determined by the physical energet- 
ics of the system; the maximum kinetic energy as de- 
termined by the energetics of the constructed surface is 
slightly higher. 

In Fig. [SI we plot the cross section as a function of H~ 
kinetic energy for several values of incident electron en- 
ergy. To compare with the experimental results of Belie, 
Landau and Hall Q, which reflect the finite resolution 
of the kinetic energy of the anion recoil, we incorporate 
the experimental resolution of 150meV in the ion kinetic 
energy direction. Such resolution effectively smears each 
vibrational peak into the next, and there are no hard ze- 
roes visible in the data of Fig. [(5] A key result of our 
calculation is that with better experimental resolution, 
the individual vibrational peaks should be able to be re- 
solved, not only for H 2 0, but also for D 2 0, and that these 
peaks should be fully separated. The experimental reso- 
lution of Ref. was insufficient to delineate the separate 



vibrational peaks for D2O. We doubt that these authors 
have resolved the rotational structure for H~+OH pro- 
duction at 7.5eV incident electron energy, as they claim. 

The isotope effect observed for this channel has 
been a matter of some interest. Compton and 
Christophoreau || observed the ratio of peak heights 
for the deuterated (D2O) to the nondeutcratcd (H2O) 
species to be 0.75, and the ratio of the energy-integrated 
cross sections, which approximate the ratio of survival 
probabilities P surv calculated with Eq. ([55)) . to be 0.60. 
In contrast, we observe a larger peak for the deuterated 
species, and similar survival probabilities P surv , both 
near 1. 

The recent results of Fedor et al. [2(| resolve this 
discrepancy. The peak heights which they obtain for 
H~+OH production versus D~+OD production via the 
2 Bi resonance indicate a larger peak for D~+OD, re- 
versing the prior experimental evidence, and putting ex- 
periment and theory on qualitatively similar ground. It 
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is clear that the ratio of peak heights obtained by Fedor 
et at, while not explicitly calculated by these authors, is 
nearer to 1 than the present theoretical results, but it is 
reassuring that the trend for both experiment and theory 
is in the same direction. The combination of the results 
of Ref. [13] and the present results indicate that the sur- 
vival probability for the physical 2 B\ state is indeed near 
1, and that minimal flux is lost via the autodetachment 
for DEA via this resonance. 

Plots of the prop agated wave packet are shown in the 
EPAPS archive (gflj- 



2. Production of H2+O via the B 1 state 

This channel is by far the minor channel for DEA via 
the 2 B\ resonance. The peak of the H2+O - cross sec- 
tion is approximately l/40 th the height of the peak for 
the major H-+OH channel @. Being such a minor chan- 
nel, it presents a more difficult challenge for theoretical 
methods such as MCTDH, and a greater test for the lo- 
cal complex potential model. We were able to obtain 
converged cross sections with the present treatment, al- 
though the magnitudes of our calculated values are far 
below the experimental results. Therefore, it is clear that 
we have not represented the dynamics into this channel 
accurately. It is possible that small discrepancies in our 
calculated surface are to blame, or that the LCP model 
is inadequate. 

We present the cross sections calculated for H2+0~ 
production as a function of incident electron energy in 
Fig. [7J We compare the total cross section for H 2 pro- 
duction from H2O with that for D2 production from D2O 
in the top panel of this figure. The cross sections peak 
at 7.6cV, 0.5eV above the experimental peak at 7.1eV, 
and are far smaller than the experimental result. Al- 
though our representation of the nuclear dynamics lead- 
ing to this channel is clearly lacking, we performed addi- 
tional calculations in which the target state of H2O was 
rovibrationally excited. We performed two calculations 
for total angular momentum J = 5, employing the cen- 
trifugal sudden (CS) approximation with K = 0, in the 
i?-embedding coordinate system, as well as a calculation 
with J = but one quantum of bend, the (010) state. 
The total cross sections for production of H2 from these 
excited states are compared to the ground initial state 
result in the top panel of Fig. [7] 

As is clear from these results, initial excitation of the 
target may play a large role in determining the magni- 
tude of the DEA cross section for H2+O - production via 
the 2 Bi state, but is insufficient to explain the discrep- 
ancy between the theoretical and experimental results. 
The effect of bending excitation increases the cross sec- 
tion dramatically; rotational excitation to J = 5 also 
enhances the cross section by approximately a factor of 
2. The excitation energy of the bending mode is approx- 
imately 0.2eV; that of the J = 5, K = state is ap- 
proximately 0.056cV. These quantities may be compared 
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FIG. 7: Cross sections calculated for production of H2+O - 
via the 2 Bi state as a function of incident electron energy. 
Top, calculated isotope effect: comparison of H2O and D2O. 
Bottom, effect of target excitation: ground initial state re- 
sult is compared to result from (010) target state with one 
quantum of bend and to result from J = 5, K = 0. 



to the value of kT at 373.15° K, which is 0.032cV. This 
comparison indicates that the bending state is not sig- 
nificantly populated in typical experimental setups and 
should not be responsible for the magnitude of the ob- 
served cross sections. Comparison of the rotational en- 
ergy to kT indicates that the degree of rotational exci- 
tation of the target may determine the precise value of 
the peak cross section observed in experiment. However, 
rotational excitation of the sample is insufficient to ex- 
plain the discrepancy between our results and the exper- 
imental ones. For the i?-embcdding coordinate system, 
the K = projection of angular momentum is the most 
likely to enhance the DEA cross section for production 
of H2, because that projection minimizes the centrifugal 
potential in the thh coordinate [see Eq. (|T7|) ] . The calcu- 
lated enhancement is due to the effect of the centrifugal 
potential in R, which "pushes" the wave packet toward 
large R, where the _ +H2 potential well lies. 

Because the width T of the 2 B\ resonance is small for 
all nuclear geometries, one might expect any nonlocal 
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effects in the resonant nuclear dynamics to be small as 
well. However, we are here considering a minor channel, 
which is only barely accessible with LCP dynamics on the 
constructed potential-energy surface. If nonlocal effects 
were to open a new dynamical pathway, or otherwise ef- 
fectively lower the dynamical barrier to the H2+0~ well, 
the magnitude of such effects would not have to be great 
in order to produce a noticeable enhancement of such a 
small cross section. Therefore, we regard nonlocal effects 
to be a strong candidate for the source of the experimen- 
tally observed cross section for production of H 2 +0~ via 
the 2 Bi resonance. 



B. Dissociative electron attachment via the 2 Ai (1 
2 A') state, Renner-Teller coupled to the 2 Bi state 

These calculations arc performed in the diabatic (l z i = 
±1) basis, which diagonalizes the nuclear kinetic energy 
operator with the Renner-Teller effect, as per the discus- 
sion in Sec. IIV( employing the centrifugal sudden Hamil- 
tonian of Eq. (fT7| . The initial state is the adiabatic 2 A\ 
(1 2 A') state, comprised of equal parts l z > = ±1. Like the 
other calculations which we present that incorporate ro- 
tational motion, they are parametrized by the body-hxed 
angular momentum quantum number K, which is the 
projection of the total angular momentum onto the em- 
bedding axis. However, for these Renner-Teller calcula- 
tions K is interpreted as the eigenvalue of the projection 
of the nuclear rotational angular momentum R z i , not the 
total angular momentum J z i, upon the embedding axis, 
and therefore the diabatic basis l z i = ±1 corresponds to 
J z , =K±1, [6l| 

We have obtained cross sections for the major H _ +OH 
(X 2 n) channel of DEA via the 2 Ai resonance. How- 
ever, for the minor H2+O - channel, we have not been 
able to obtain converged, nonzero cross sections. The 
mechanism for DEA via the 2 A x resonance to produce 
H2+O - remains unknown. It is possible that three-body 
breakup, which we have not treated, is important here. 

The considerations of Rcf. [l8[ indicate that the nu- 
clear dynamics of DEA via the 2 Ai resonance may hold 
some surprises, and that the LCP model may be insuf- 
ficient for a full description thereof. In particular, as 
demonstrated in paper I, the width of the 2 A\ resonance 
becomes large as the nuclear geometry moves toward the 
H~+OH product arrangment. We have calculated width 
values as high as 0.15cV for this resonance state for such 
stretched geometries, despite the fact that the resonance 
state lies only slightly above the neutral at these geome- 
tries and ultimately becomes bound as the atom-diatom 
distance is further increased. The explanation for this be- 
havior is that the electronic structures of the neutral and 
anion become highly correlated and different from each 
other at such stretched geometries, and as a result, there 
is considerable shape resonance character mixed into the 
2 A\ Feshbach resonance. The radically peaked behavior 
of the width of the 2 A\ state may portend a breakdown of 



the LCP model, which relics on the implicit assumption 
that the background-resonance coupling is a relatively 
smooth function of nuclear geometry. Also, the fact that 
the 2 A\ resonance is coupled to the neutral target by an 
s-wave matrix element indicates that virtual state effects 
may play a large role as the resonance becomes bound. 
Such virtual state effects cannot be properly described 
by the LCP model, but have been taken into account in 
other systems using effective range theory as, for exam- 
ple, in Refs. [13, Hj|. The fact that we have overesti- 
mated the magnitude of the experimental cross section 
for H~+OH production via the 2 A\ resonance indicates 
that a breakdown of the LCP model may be responsible 
for the loss of flux via autodctachment. 



1. Production of H + OH via the Renner-Teller coupled 
2 A" CBi) & 1 2 A' ( 2 A 1 ) states 

The cross sections for total H _ +OH production were 
easy to converge, as the channel involved is the main 
channel, the dissociation direct, and the dynamics ap- 
parently reasonably separable in the toh Jacobi internal 
coordinate system. The relatively small size of the single 
particle function (SPF) expansions required to converge 
the calculation (see the EPAPS archive[6(|) support this 
conclusion. We found the Renner-Teller coupling to have 
a negligible effect on both the magnitude of the total cross 
section and its breakdown into rotational and vibrational 
states, so we only report results for K = 0. The cross sec- 
tions (for K — 0) into final vibrational and rotational lev- 
els of the diatomic fragment, for both the deuterated and 
nondeuterated cases, are given in the EPAPS archive [60|. 

Our theoretical treatment overestimates the cross sec- 
tion for DEA into this channel via the 2 A\ resonance. 
Our peak heights, 4.14 and 4.16 xlO" 18 cm 2 for the 
nondeuterated and deuterated target, respectively, are 
approximately three times larger than Melton's observed 
peak height of 1.3 xlO -18 cm 2 . Although the recent re- 
sults of Fedor et al. [Icj do resolve this peak better than 
previous experiments, and indicate that Melton's peak 
height value may be too low, there is a clear discrep- 
ancy between theory and experiment here. We attribute 
the discrepancy to virtual state effects, as discussed in 
Ref. [lH, which may lead to significant autodetachment 
as the 2 Ai (1 2 A') state becomes bound. There is also 
the possibility that we have overestimated the entrance 
amplitude for this state, as we have done for the 2 B\ 
state. 

The degree of vibrational excitation is higher for this 
resonance than for the B\ resonance: the values of (u) 
calculated from Eq.® are 2.03 and 1.28 for the 2 A x 
and 2 B\ resonances, respectively. This difference is most 
likely due to the gradient of the potential-energy surface 
in the symmetric stretch direction, which is larger at the 
equilibrium geometry of the neutral for the 2 A\ surface 
than for the 2 B\ surface. The behavior of the propagated 
wave packet, which is plotted in the EPAPS archive [60|. 
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FIG. 8: Top, production of H~+OH from the 2 Ai state, as 
a function of incident electron energy and H~ fragment ki- 
netic energy, as in previous plots; bottom, deuterated. The 
maximum kinetic energy available, as determined from the 
physical energetics, is plotted with a bold line. Contours ev- 
ery 6 x 10" 18 cm 2 eV -1 . 
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FIG. 9: Production of H"+OH ( 2 I1) and D~+OD ( 2 I1) via 
the 2 Ai state at 8.5eV incident electron energy, as a function 
of fragment kinetic energy, as in previous plots. Calculated 
results have been broadened using a 150meV linewidth. to 
compare with the experimental results from Belie, Landau 
and Hall also plotted. 



final anion recoil kinetic energy. Fig. [H] displays a clear 
difference from Fig. This figure shows that for DEA 
via the 2 A\ state, the degree of rotational excitation of 
the diatomic fragment is high enough that the different 
vibrational states are distinguishable, but not completely 
separated. 

We compare our results for the laboratory-frame, anion 
recoil kinetic energy distribution at an incident electron 
energy of 8.5eV with the corresponding results of Belie, 
Landau and Hall@, in Fig. [5] (Two-dimensional plots of 
the cross section as a function of both incident electron 
energy and final anion recoil kinetic energy can be found 
in the EPAPS archive [60l].) In contrast to the case of H~ 



is similar to that found for the 2 B\ resonance: the wave 
packet experiences an initial impulse in the symmetric 
stretch direction, but then is bifurcated by the developing 
potential wall in this direction, and reflected into cither 
H~+OH channel. The vibrational excitation is the result 
of the wave packet oscillating in the r direction as it 
passes down the OH potential well. 

The degree of rotational excitation within the OH frag- 
ment is also higher for the 2 A\ state than for the 2 B\ 
state. Using Eq. (pS5)) . we calculate an average degree of 
rotational excitation <j 2 ) = 119 for this resonance, com- 
pared to 107 for the 2 B\ resonance. This results from 
the larger gradient of the potential-energy surface in the 
bend direction for the 2 A\ surface compared to the 2 B\ 
surface. The 2 A\ wave packet is given an impulse in the 
bend direction, which corresponds to excitation of rota- 
tional quanta j. This excitation persists in the final state, 
as demonstrated by these calculations. 

Fig. [5] presents two-dimensional plots of the cross sec- 
tion as a function of both incident electron energy and 
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FIG. 10: Production of D"+OD ( 2 II) at 9.3eV ( 2 Ai) and 
6.8eV ( 2 Bi), as a function of D fragment kinetic energy, as 
in previous figures, with experimental results from Curtis and 
Walker []7|. 
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production from the 2 Bi state, here the degree of rota- 
tional excitation of the diatomic fragment is high enough 
that the different vibrational states are distinguishable, 
but not completely separated. The experimental resolu- 
tion of Ref. Q was insufficient to delineate the different 
peaks for various OH (i/), although Fig. [8] demonstrates 
that with sufficient resolution, the vibrational structure 
should be apparent. In Fig. 1101 we compare to the ex- 
perimental results of Curtis and Walker [lOj. 



2. Failure to calculate production of H2 + O via 
dissociative attachment to the 2 A\ state 

We have been unable to obtain a nonzero cross sec- 
tion for DEA via the 2 A\ (1 2 A') state, Renner- Teller 
coupled to the 2 B\ state, leading to H 2 +0 _ . We have 
attempted calculations for K=0 (uncoupled), 1, 2, 3, and 
4. Within the MCTDH calculations, we employed single- 
particle function (SPF) expansions of up to 24x36x30, 
with no success. With this large SPF expansion, and 
propagation times of up to 100 fs, we regard the rep- 
resentation of the LCP dynamics within the MCTDH 
ansatz to be accurate. We suspect that O - production 
from 2 Ai may be dominated by three-body breakup into 
H+H+0~, which we have not treated. 



C. Dissociative electron attachment via the 2 £>2 
state, involving the conical intersection with the 2 Ai 

state 

As described in Ref. [H[ and paper I, dissociative elec- 
tron attachment to H2O via the highest-energy 2 B^ state 
must involve the conical intersection that this state ex- 
hibits with the 2 A\ state. The gradient of the potential- 
energy surface leads directly toward this conical intersec- 
tion from the equilibrium geometry of the neutral. The 
conical intersection forms a line in the three-dimensional 
space of nuclear geometries, and occurs within C2V sym- 
metry, where the OH bond lengths are equal. 

We performed a diabatization on the results of 
configuration-interaction calculations on 2 A\ and 2 Bi (1 
and 2 2 A') resonances, as described in paper I, to pro- 
duce a set of diabatic 2 A\ and 2 Bi surfaces along with 
a coupling surface. These diabatic surfaces are employed 
in all of the following calculations. 

Before describing the individual calculations, a few 
preliminary remarks about the experimental observations 
are in order. Although absolute cross sections for H~ 
production via the 2 i?2 resonance are not available, the 
experimental evidence [H,[2(j indicates that for both D2O 
and H2O target states, the branching ratio between H~ 
production and 0~ production highly favors 0~. There- 
fore, the dynamics of DEA beginning in the 2 B^ state are 
much different from those for the lower-energy 2 A\ and 
2 B\ resonances, which yield far more H _ than O - . 

This observation is not surprising, in light of the 



potential-energy surfaces which we have calculated and 
shown in paper I; the upper 2 2 A! surface was demon- 
strated to be quite different from those of the lower reso- 
nances. In particular, the dynamics beginning on the 2 Bi 
(2 2 A') surface will begin with a decrease of the H-O-H 
bond angle Ohoh, which motion will favor the H2+0 - 
product arrangement. However, as we will show, there 
appears to be active competition between the two prod- 
uct arrangements, and the branching ratio observed in 
experiment is likely the product of both the shape of the 
real- valued component to the 2 2 A' surface and the be- 
havior of its imaginary component. 

As is clear from Table HI the cross sections we calcu- 
late for this channel are smaller than the observed cross 
sections. However, the comparision is complicated by 
the fact that the three-body dissociation channels are 
open for incident electron energies sufficient to reach the 
2 i?2 resonance, and we produce cross sections only for 
the two-body dissociation channel; the disagreement may 
therefore be due to a large contribution of the three- 
body breakup channel to the dominant production of 
0~. However, the locations of the calculated and ex- 
perimental peak maxima for production of O - from the 
2 £?2 resonance agree very well: both cross sections peak 
at about 11.8cV. Although the presence of the three- 
body dissociation channel may shift the peak, this com- 
parison indicates that we have probably accurately rep- 
resented the vertical transition energy for the 2 Bi res- 
onance. The vertical transition energy as represented 
by our configuration-interaction surface is 12.83eV, and 
therefore we recommend a value of approximately 12.8eV 
for the appropriate physical transition energy. This value 
is above the value of 11.97cV given by the complex Kohn 
calculations of paper I. 

The calculated branching ratio between H2+0 - 
(D 2 +0~) and OH+H~ (OD+D - ) production is near 
unity, but the experimental ratio (for the undeuterated 
product) exceeds 1 by a large factor. As we will show, the 
dynamics within the first few femtoseconds after attach- 
ment are controlled by both real and imaginary compo- 
nents of the potential-energy surface, the latter consum- 
ing most of the propagated wave packet within the first 
six femtosectonds. 

We have examined the effect of rotational excitation 
upon the cross section for production of H2+O - from 
DEA via the 2 £>2 state, and find it to be negligible. 



1. Production of H~ + OH and 2 T,) via the upper B2 

state 

The calculation for the production of OH ( 2 n and 
2 E)+H~ via the 2 B^ Feshbach resonance, which is cou- 
pled to the 2 A\ resonance via their conical intersection, 
proved difficult to converge. This is evidenced by ragged- 
ness in the OH ( 2 n) channel cross sections. A final state 
resolution in this channel was not possible , although we 
were able to resolve the final states of OH ( 2 S). 
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FIG. 11: Top, total cross section calculated for production of 
H~+OH ( 2 n) versus ET+OH ( 2 E) from 2 B 2 state as a func- 
tion of incident electron energy. Bottom, deuterated version. 



With reference to the discussion in Ref. [l8( and paper 
I, the constructed diabatic 2 B 2 surface correlates to the 
species H~+OH ( 2 E), as does the adiabatic 2 B 2 (2 2 A') 
surface. The diabatic 2 A\ surface correlates to H~ plus 
ground state OH ( 2 II), as does the adiabatic 1 2 A' state. 
Therefore, dynamics beginning on the 2 B 2 surface that 
leads to production of the ground-state H~+OH ( 2 IT) 
species must proceed via the off-diagonal coupling to the 
2 A\ surface. From the viewpoint of the diabatic basis, 
the off-diagonal coupling must in this case lead to a tran- 
sition between the diabatic 2 B 2 and 2 A\ surfaces; from 
the viewpoint of the adiabatic basis, the dynamics must 
proceed through the conical intersection via the singular 
derivative couplings inherent in that basis. 

We present the calculated total cross sections for pro- 
duction of either H~+OH or D~+OD in Fig. [TTJ The 
results are similar in shape, but the magnitude of the 
cross sections for the deuterated case are approximately 
half those of the nondeuterated case. Differences in 
the reduced masses in the dissociative direction result 
in a relatively longer time during which the deuterated 
species may undergo autodetachmcnt and, consequently, 
a smaller survival probability for the deuterated anion, 



As shown in TablcHH t ne survival probability for the non- 
deuterated 2 B 2 state is 21.5%, whereas for the deuterated 
species it is only 13.1%. Unfortunately, there is no ex- 
perimental data for comparison that measures the rela- 
tive magnitude of the H~ and D~ peaks for the highest- 
energy 2 B 2 resonance. 

An obvious feature of the results presented in Fig. [TTlis 
that the branching ratio of OH ( 2 n) to OH ( 2 S) produc- 
tion depends on the incident energy of the electron. This 
ratio varies from 100% (only 2 n is produced) at onset to 
zero (only 2 S produced) at higher energy. At low energy 
the observed cross section is the result of dynamics in 
which the wave packet makes a nonadiabatic transition 
from the upper 2 B 2 (2 2 A') surface to the lower 2 A\ (1 
2 A') surface, whereas at high energy, the observed cross 
section is due to dynamics in which there is no transition. 
Thus, the nuclear dynamics via the 2 B 2 state involve the 
conical intersection to produce a branching ratio that 
varies with incident energy in an interesting way. 

We have only been able to achieve final-state resolu- 
tion for the OH ( 2 E) fragment. Tow-dimensional views 
of the cross sections for H~+OH( 2 £) production, as a 
function of both the incident electron energy and the ki- 
netic e nerg y of the H~ fragment, are given in the EPAPS 
archive 60], along with comparisons of our calculated re- 
sults with previous experiment. Our calculations repro- 
duce the approximate level of excitation within the di- 
atomic fragment, as the theoretical and experimental re- 
sults are both centered near the same kinetic energy, 
~2.75eV for H~ from H2O, and ~1.5eV for D~ from 
D2O. We cannot make a more quantitative comparison, 
because there are no experimental values for the average 
kinetic energy release in this channel. 

The wave packet dynamics for DEA leading to H~ pro- 
duction via the 2 B 2 state coupled to the 2 A\ state are 
shown in Figs. [T2l and IT51 The former shows the reduced 
density on the adiabatic 2 2 A' surface, where it starts; the 
latter shows that on the adiabatic 1 2 A' surface. These 
plots were obtained by transforming the propagated wave 
packets from the diabatic basis to the adiabatic basis. 
The wave packet initially has no magnitude on the lower 
1 2 A' surface. Nonadiabatic coupling changes this situa- 
tion as the wave packet is propagated. The norm of the 
propagated wave packet on the 1 2 A' surface reaches a 
maximum of 0.112 at £=9.1 fs by which time a portion 
of the wave packet has reached the dissociative H _ +OH 
( 2 n) well of the 1 2 A' surface. The portion of the wave 
packet within this well (see the bottom-left panel of Fig. 
[T3")) lies beyond i?=4.5an where the resonance becomes 
bound, and so it continues toward dissociation with neg- 
ligible loss of flux. The subsequent decrease of the norm 
of the wave packet on the 1 2 A' surface is therefore due 
to the consumption of other parts of the wave packet 
by the imaginary component of this surface, and to its 
absorption by the complex absorbing potentials. 

As described in paper I, the magnitude of the width 
for the upper 2 2 A' surface is generally large, though it 
decreases slowly as the H~+OH ( 2 E) well is approached, 
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FIG. 13: Propagation of wave packet on coupled 2 £?2 and 2 A\ surfaces for H~+OH ( 2 II / 2 E) channels, adiabatic 1 2 A' (— > 
2 II) component, with real part of 1 2 A' potential-energy surface at 7=90°. Bond lengths, units of bohr. Density is integrated 
over 7. 



and abruptly as the Eb+O - well is approached. As a 
result, the wave packet which begins upon the upper 2 
2 A' surface is rapidly consumed, and its norm decreases 
from exactly 1 to 0.321 within six femtoseconds. At this 
time, the combined norm on both surfaces is 0.402. The 
calculated total survival probability for this resonance, 



P S urv, calculated with Eq. ([3lj) . is 21.5% (see Table [n}. 
From this comparison we can see that the majority of 
the autodetachment for this resonance occurs within the 
first six femptoseconds; its survival probability is 40.2% 
within this initial time period, and 21.5/40.2 = 53.5% 
thereafter. 
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FIG. 14: Cross sections for production of H2 (v)+0~ from 
2 E>2 state, with rotational excitation of the target, as a func- 
tion of incident electron energy. Bold line, total cross sections 
for ground rotational state (J = 0); thin line, calculations for 
J = 5, K = 0. 

The high degree of vibrational excitation ((v)) in both 
the OH (II) and OH (£) channels is apparent in the os- 
cillations of the dissociating wave packet within each po- 
tential well, visible in the lower panels of Figs . fT2l and [T3l 
Fig. [13] shows that in this case there is additional struc- 
ture to the dissociating wave packet on the lower 1 2 A' 
surface; however, this structure is most likely due to the 
calculation not being fully converged. 

2. Production of H2 + 0~ via the upper 2 i?2 state 

The channel H 2 +0 _ is the dominant channel observed 
in experiment for dissociative attachment to water via 
the highest-energy 2 B 2 resonance. As discussed at length 
in Ref. (l8j and paper I, this channel is not present as an 
asymptote on the 2 B 2 (2 2 A') surface, and therefore, the 
system must undergo a nonadiabatic transition via the 
conical intersection to the lower 1 2 A' surface in order 
to reach this product channel. In the context of the rep- 
resentation which we constructed in paper I, the system 
must follow the diabatic 2 B 2 surface past its crossing with 
the 2 Ai diabatic surface. The H2+0~ channel is present 
as an as ymp tote of the diabatic 2 B 2 surface. As described 
in Ref. [18| and paper I, the adiabatic 2 A\ surface does 
not have a bound asymptote in this arrangement; it cor- 
relates to 0~+H2 (CgO"i) instead. 

The calculated peak cross section for this channel, 
1.87xl0 _19 cm 2 , is smaller than the experimental value, 
5.7x 10 _19 cm 2 , reported by Melton Q . The comparison 
with experiment is again complicated by the fact that we 
calculate only the two-body DEA cross section, while the 
available experimental data do not discriminate between 
production of 0~+H+H and _ +H 2 . A possible expla- 
nation for the discrepancy between our calculated results 
and experiment is the presence of a large three-body corn- 
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FIG. 15: Cross sections for production of H2 (u)+0~ from 
2 i?2 state as a function of incident electron energy, showing 
degree of vibrational and rotational excitation of the H2 frag- 
ment. Total cross section, thick line. Top, cross sections 
summed over vibrational quantum number v and binned by 
rotational quantum number j; bottom, v and j reversed. 



poncnt to 0~ production via the 2 B 2 resonance. Rota- 
tional excitation of the target H2O molecule, on the other 
hand, cannot account for this discrepancy. We have per- 
formed several calculations in which rotational excitation 
of the target is included. These include calculations for 
total angular momentum J — 5, projection K = 0. We 
find that the effect of such rotational excitation is mini- 
mal, as Fig. [T4l shows. 

We calculate a very high degree of rotational and vi- 
brational excitation in the H2 or D 2 fragment. The aver- 
age degree of vibrational excitation (y) calculated from 
Eq.([3S]is 7.75 for the H 2 fragment and 13.0 for the D 2 
fragment. The corresponding values for (j 2 ) are 405 and 
725, respectively. Figure [15] shows the total cross sec- 
tions, as well as the cross sections into either rotational 
or vibrational states, summed over the opposite quantum 
number. The degree of vibrational excitation evidently 
decreases with incident electron energy, while the degree 
of rotational excitation shows little correlation with inci- 
dent electron energy. 
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FIG. 16: Propagation on coupled 2 i?2 and 2 Ai surfaces for Hh+O - channel, adiabatic 2 2 A' component. The reduced density 
(integrated over 7) of the adiabatic 2 2 A' component of the propagated wave packet is plotted with the real part of the 2 2 A' 
potential-energy surface at 7 = 90° (C21) geometry). The location of the conical intersection is marked with a bold line. Bond 
lengths, units of bohr. 



The high degree of rotational and vibrational excita- 
tion of the diatomic fragment reduces the kinetic energy 
of the atom-diatom recoil. This is reflected in the cross 
sections for production of both H2 and D2 via the 2 i?2 
resonance, which have the greatest magnitude nearere 
the lower range of recoil energy. Two-dimensional plots 
of these cross sections as functions of both incident elec- 
tron energy and the kinetic energy of the recoil are shown 
in the EPAPS archive^. 

Plots of the propagated wave packet for DEA via the 
2 £?2 resonance are shown in Figs. [TBI and [T71 The first of 
these shows the magnitude-squared of the 1 2 A' compo- 
nent to the propagated wave packet, integrated over 7, 
and the latter shows the 2 2 A' component. The corre- 
sponding potential-energy surfaces, evaluated at 7=90°, 
are also plotted, along with the location of the conical 
intersection seam which appears as a bold line. 

The wave packet begins on the upper surface and pro- 
ceeds to the lower surface only via nonadiabatic coupling 
near the conical intersection. As described in paper I, 
the gradient of the upper 2 2 A' resonance surface leads 
downhill toward its conical intersection with the 1 2 A' 
resonance, leading the propagated wave packet toward 
the seam. This behavior is clearly visible in Fig. [16] and 
[T71 The 2 2 A' wave packet follows the conical intersection 



seam in Fig. II 61 until it is consumed by the large imagi- 
nary component to that potential-energy surface and by 
nonadiabatic coupling to the 1 2 A' state along the inter- 
section. The wave packet appears on the 1 2 A' surface in 
FigfTT] along the conical intersection, and a small portion 
of it is able to reach the H2+O - well of that state. 



The magnitude of the cross section for production of 
H2+O - from the 2 Bi resonance is therefore controlled by 
several competing effects. The shape of the real part of 
the potential-energy surface determines the dynamically 
accessible pathways and favors localization of the 2 2 A' 
wave packet near the conical intersection. At the same 
time, the large imaginary component to this surface con- 
sumes the wave packet and decreases the amount of flux 
available to enter the conical intersection. On the lower 
1 2 A' surface, the amount of flux that enters the H2 po- 
tential well is determined by the shape of that potential- 
energy surface, since the conical intersection is outside 
the potential well and only a fraction of the wave packet 
is propagated into the well. 
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FIG. 17: Propagation on coupled 2 i?2 and 2 A\ surfaces for H2+0~ channel, adiabatic 1 2 A' component. The reduced density 
(integrated over 7) of the adiabatic 1 2 A' component of the propagated wave packet is plotted with the real part of the 1 2 A' 
potential-energy surface at 7=90° (Civ geometry). The location of the conical intersection is marked with a bold line. 



VII. CONCLUSION 

We have presented the results of a fully ab initio study 
of dissociative electron attachment to H2O that includes 
the full dimensionality of nuclear motion. We have at- 
tempted to calculate the cross sections for all the ma- 
jor and minor two-body channels which are present as 
asymptotes of the Born-Oppcnheimer, 2 Bi, 2 Ai, and 
2 i?2 adiabatic electronic Feshbach resonances. While we 
have qualitatively described the principal features that 
have been experimentally observed, it is clear that a 
fully quantitative description of this process has yet to 
be achieved. 

The nuclear dynamics calculations were carried out us- 
ing the MCTDH method within the framework of the 
local complex potential model. For the major channel 
DEA, H^+OH (X 2 II) production through the lowest 
2 Bi resonance, the underlying assumptions of the model 
are well satisfied and we have obtained reasonably good 
agreement with the experimental observations. Another 
notable feature of the present study is the quantification 
of the mechanism in the major channel that leads to pro- 
duction of O - through the 2 Bi resonance. Our earlier 
speculation [l8| that a conical intersection between the 
2 A\ and 2 £?2 states would play the key role in this process 
has been confirmed by the present study. 



The present treatment has been limited to a consid- 
eration of DEA only into the final-state two-body chan- 
nels. This limitation undoubtedly explains our inability 
to produce a non-zero cross section for 0~ production 
via the 2 A\ resonance, which is likely to be dominated 
by three-body breakup. Three-body breakup may also 
play a role in 0~ production via the 2 i?2 resonance, and 
its neglect here could explain why our calculated cross 
sections are smaller than the experimental results, which 
did not differentiate two- and three-body channels. 

Physics beyond the local complex potential model may 
be at work in some of the minor channels. Dissocia- 
tive electron attachment via the 2 A\ Feshbach resonance 
may involve an even greater variety of complicated res- 
onant as well as non-resonant phenomena not described 
by the LCP model. A variety of effects that go beyond 
the LCP model could be at play in the production of 
H~ via 2 A\ Feshbach resonance, including coupling to 
a broader shape resonance and even non-resonant vir- 
tual state effects. The neglect of such effects could well 
explain our ovcrcstimation of the cross section for pro- 
duction of H~+OH via the 2 A\ state. Even for DEA 
via the lowest-energy 2 B\ state, nonlocal physics may be 
important in the minor channel, which leads to H2+0~. 

We have achieved considerable success in describing 
the mean features of DEA to water, clarified the mech- 
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anisms for the two-body breakup channels, and found 
evidence to suggest that three-body breakup to produce 
0~ might be important. Nonetheless, many challenges 
remain before a complete and quantitative understand- 
ing of this fundamental, but complicated, system will be 
realized. 



nia Lawrence Berkeley National Laboratory under Con- 
tract DE-AC02-05CH11231 and was supported by the 
U.S. DOE Office of Basic Energy Sciences, Division of 
Chemical Sciences. The authors acknowledge many help- 
ful discussions with H. -Dieter Meyer concerning various 
aspects of the MCTDH method. 



Acknowledgments 



This work was performed under the auspices of the 
US Department of Energy by the University of Califor- 



[1] D. J. Haxton, C. W. McCurdy, and T. N. Rescigno, Phys. 

Rev. A 75, 012711 (2007). 
[2] W. N. Lozier, Phys. Rev. 36, 1417 (1930). 
[3] I. S. Buchel'nikova, Zh. Eksperim. i Teor. Fiz. 35, 1119 

(1959). 

[4] G. J. Schultz, J. Chem. Phys. 44, 3856 (1966). 

[5] R. N. Compton and L. G. Christophorou, Phys. Rev. 

154, 110 (1967). 
[6] C. E. Melton, J. Chem. Phys. 57, 4218 (1972). 
[7] L. Sanche and G. J. Schultz, J. Chem. Phys. 58, 479 

(1972). 

[8] S. Trajmar and R. I. Hall, J. Phys. B. 7, L458 (1974). 
[9] D. S. Belie, M. Landau, and R. I. Hall, J. Phys. B. 14, 
175 (1981). 

[10] M. G. Curtis and I. C. Walker, J. Chem. Soc. Faraday 

Trans. 88, 2805 (1992). 
[11] C. R. Claydon, G. A. Segal, and H. S. Taylor, J. Chem. 

Phys 54, 3799 (1971). 
[12] M. Jungen, J. Vogt, and V. Staemmler, Chem. Phys. 37, 

49 (1979). 

[13] T. J. Gil, T. N. Rescigno, C. W. McCurdy, and B. H. 

Lengsfeld HI, Phys. Rev. A 49, 2642 (1994). 
[14] L. A. Morgan, J. Phys. B 31, 5003 (1998). 
[15] J. D. Gorfinkel, L. A. Morgan, and J. Tennyson, J. Phys. 

B 35, 543 (2002). 
[16] D. J. Haxton, Z. Zhang, C. W. McCurdy, and T. N. 

Rescigno, Phys. Rev. A 69, 062713 (2003). 
[17] D. J. Haxton, Z. Zhang, H.-D. Meyer, T. N. Rescigno, 

and C. W. McCurdy, Phys. Rev. A 69, 062714 (2003). 
[18] D. J. Haxton, T. N. Rescigno, and C. W. McCurdy, Phys. 

Rev. A 72, 022705 (2005). 
[19] D. J. Haxton, T. N. Rescigno, and C. W. McCurdy, Phys. 

Rev. A 73, 062724 (2006). 
[20] J. Fedor, P. Cicman, B. Coupier, S. Feil, M. Winkler, 

K. Gluch, J. Husarik, D. Jaksch, B. Farizon, N. J. Mason, 

et al., J. Phys B 39, 3935 (2006). 
[21] D. T. Birtwistle and A. Herzenberg, J. Phys. B 4, 53 

(1971). 

[22] L. Dube and A. Herzenberg, Phys. Rev. A 20, 194 (1979). 
[23] J. N. Bardsley and J. M. Wadehra, J. Chem. Phys. 78, 
7227 (1983). 

[24] T. F. O'Malley and H. S. Taylor, Phys. Rev. 176, 207 
(1968). 

[25] T. F. O'Malley, Phys. Rev. 150, 14 (1966). 
[26] C. E. Klots and R. N. Compton, J. Chem. Phys. 69, 1644 
(1978). 



[27] S. Goursaud, M. Sizun, and F. Fiquet-Fayard, J. Chem. 

Phys 65, 5453 (1976). 
[28] S. Goursaud, M. Sizun, and F. Fiquet-Fayard, J. Chem. 

Phys 68, 4310 (1978). 
[29] H. Feshbach, Ann. Phys. 19, 287 (1962). 
[30] A. U. Hazi, T. Rescigno, and M. Kurilla, Phys. Rev. A 

23, 1089 (1981). 
[31] J. Bieniek, J. Phys. B 13, 4405 (1980). 
[32] C. W. McCurdy and J. L. Turner, J. Chem. Phys. 78, 

6773 (1983). 

[33] J. Tennyson and B. T. Sutcliffe, J. Chem. Phys. 77, 4061 
(1982). 

[34] J. Z. H. Zhang, Theory and Application of Quantum 
Molecular Dynamics (World Scientific Publishing Co. 
Pte. Ltd., 1999). 

[35] A. R. Edmonds, Angular Momentum in Quantum Me- 
chanics (Princeton University Press, 1996). 

[36] C. Petrongolo, J. Chem. Phys. 89, 1297 (1988). 

[37] S. Sukiasyan and H.-D. Meyer, J. Phys. Chem. A 105, 
2604 (2001). 

[38] T. Mulloney and G. C. Schatz, Chem. Phys. 45, 213 
(1980). 

[39] R. T. Pack, Chem. Phys. Lett. 108, 333 (1984). 

[40] R. Renner, Z. Phys. 92, 172 (1934). 

[41] C. Jungen and A. J. Merer, Mol. Phys. 40, 1 (1980). 

[42] S. Carter and N. C. Handy, Mol Phys. 52, 1367 (1984). 

[43] A. Loettgers, A. Untch, H.-M. Keller, R. Schinke, H.-J. 

Werner, C. Bauer, and P. Rosmus, J. Chem. Phys. 106, 

3186 (1997). 

[44] I. B. Bersuker, Chem. Rev. 101, 1067 (2001). 

[45] A. V. Fishchuk, P. E. S. Wormer, and A. van der Avoird, 

J. Phys. Chem. A 110, 5273 (2006). 
[46] A. V. Fishchuk, G. C. Groenenboom, and A. van der 

Avoird, J. Phys. Chem. A 110, 5280 (2006). 
[47] G. C. Schatz and A. Kuppermann, J. Chem. Phys. 65, 

4642 (1975). 

[48] H.-D. Meyer, U. Manthe, and L. S. Cederbaum, Chem. 

Phys. Lett. 165, 73 (1990). 
[49] U. Manthe, H.-D. Meyer, and L. S. Cederbaum, J. Chem. 

Phys. 97, 3199 (1992). 
[50] M. Beck, A. Jackie, G. Worth, and H.-D. Meyer, Phys. 

Reports 324, 1 (2000). 
[51] H.-D. Meyer and G. A. Worth, Theor. Chem. Acc. 109, 

251 (2003). 

[52] G. A. Worth, M. H. Beck, A. Jackie, and H.- 
D. Meyer, The MCTDH Package, Version 8.2, 



25 



(2000) . H.-D. Meyer, Version 8.3 ( 2002). See 



|http:/ / www.pci.uni-heidelberg.de /tc /usr/mctdh/ 
[53] J. C. Light, in Time-Dependent Quantum Molecular Dy- 
namics, edited by J. Broechkove and L. Lathouwers 
(Plenum, New York, 1992). 
[54] C. Leforestier and R. E. Wyatt, J. Chem. Phys 78 (1983). 
[55] R. Kosloff and D. Kosloff, J. Comput. Phys. 63 (1986). 
[56] A. Jackie and H.-D. Meyer, J. Chem. Phys. 105, 6778 
(1996). 

[57] J. Rost, J. Griffin, B.Friedrich, and D. Hershbach, Phys. 

Rev. Lett. 68, 1299 (1992). 
[58] G. C. Corey and D. Lemoine, J. Chem. Phys. 97, 4115 

(1992). 

[59] J. W. T. G. C. Corey and D. Lemoine, in Numerical Grid 
Methods and Their Application to Schrddinger's Equa- 



tion, edited by C. Cerjan (Kluwer Academic Publishers, 
Netherlands, 1973), pp. 1-23. 

EPAPS Document No. E-PLRAAN-75-005702 for 
subroutines that generate the real and imaginary 
components of hte constructed potential energy 
surfaces. For more information on EPAPS, see 
http:/ /www. aip.org/pubservs/cpaps. html 



[60] 



[61] C. W. McCurdy, W. A. Isaacs, H.-D. Meyer, and T. N. 

Rescigno, Phys. Rev. A 67, 042708 (2003). 
[62] W. Vamoose, C. W. McCurdy, and T. N. Rescigno, Phys. 

Rev. A 68, 052713 (2003). 
[63] W. Vanroose, Z. Zhang, C. W. McCurdy, and T. N. 

Rescigno, Phys. Rev. Lett. 92, 053201 (2004). 



